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Abstract. We present an algorithm for the computation of the topological type of a 
real compact Riemann surface associated to an algebraic curve, i.e., its genus and the 
properties of the set of fixed points of the anti-holomorphic involution r, namely, the 
number of its connected components, and whether this set divides the surface into one 
or two connected components. This is achieved by transforming an arbitrary canonical 
homology basis to a homology basis where the ,4-cycles are invariant under the anti- 
holomorphic involution r. 



1. Introduction 

Riemann surfaces have many applications in physics and mathematics as in topological 
field theories and in the theory of integrable partial differential equations (PDEs). In 
concrete applications such as solutions of PDEs, e.g. Korteweg-de Vries and nonlinear 
Schrodinger (NLS) equations, see e.g. |3j and references therein, physical quantities as for 
instance the amplitude of a water wave are real. Thus reality conditions on the solutions 
are important in practice. The corresponding solutions have to be constructed on real 
Riemann surfaces, i.e., surfaces with an anti-holomorphic involution r acting as the complex 
conjugation on a local parameter on the surface. Regularity conditions for these solutions 
depend on the topological type of the surface, i.e., whether there are connected sets of fixed 
points of the involution r, the real ovals, and whether these ovals separate the surface into 
two connected components. 

It is well known that all compact Riemann surfaces can be realized via nonsingular 
algebraic curves in P 2 (C). F. Klein |28] observed that a real Riemann surface can be 
obtained in an analogous way from a nonsingular real plane algebraic curve 1Z with an 
affine part of the form 

N M 

(1.1) f(x, y) = J2Yl a mnX m V n = 0, x, y £ C, a mn £ R. 

n=l m=l 



The focus of this paper is on real compact Riemann surfaces. For curves of the form (1.1 ) 
the action of the complex conjugation gives rise to an antiholomorphic involution r defined 
on 1Z by r(x,y) = (x,y). The set of fixed points of r is denoted by H(M) and is called 
the real part of 1Z. The connected components of H(M) are called real ovals. Historically, 
the first result in the topology of real algebraic curves was obtained by Harnack [21]: the 



number k of real ovals of a curve 1Z of genus g satisfies < k < g + 1. In other words, the 
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number of connected components of the real part of a real nonsingular plane algebraic curve 
cannot exceed g + 1. Curves with the maximal number of real ovals are called M-curves. 

The complement 1Z\1Z(M) has either one or two connected components: if 1Z\1Z(M) has 
two components, the curve 1Z is called a dividing curve, otherwise it is called non-dividing 
(notice that an M-curve is always a dividing curve). The topological type of 1Z is usually 
denoted by (g, k, a) where g denotes the genus, k the number of real ovals, and a = if 
the curve is dividing, a = 1 if it is non-dividing. This implies that the topological type of 
a curve without real oval is (g, 0, 1). Notice that the first part of Hilbert's 16th problem is 
concerned with the relative configuration of real ovals of a plane algebraic curve of given 
degree in P 2 (M), i.e., how many ovals can lie in interior of another oval. This question 
has been studied by many authors, see, for instance, |23[ |32| I3~T| l20l [39| [T] and references 
therein. However, until now the complete answer is known only for curves of degree 7 and 
less. We will not discuss this topic here, but we would like to mention that in general, a 
solution to this problem, namely, the knowledge of the embedding 7£(K) C P 2 (M), does not 
provide any information on the embedding 7£(R) C 1Z which is the subject of the present 
paper. For instance, it is possible to construct two real plane algebraic curves having the 
same degree and the same configuration of ovals, one of them being dividing and the other 
non-dividing (see [18J p. 8). 

The aim of this paper is twofold: to determine the topological type (g, k, a) of a real 



algebraic curve of the form (1.1) with a numerical approach, and to transform periods of 
the holomorphic differentials of the curve to a form where the .A-periods are real. There 
exist various algorithms which give the oval arrangements of a given real algebraic curve, 
see, for instance, [2j El E31 HH [191 ESI EI], all of them following the same scheme. But to 
the best of our knowledge, there exists no algorithm that computes the parameter a in 
the topological type (g,k,a) of 1Z, which encodes the property of the curve to be dividing 
or not. The starting point of our algorithm is the work by Deconinck and van Hoeij 
who developed an approach to the symbolic-numerical treatment of algebraic curves. This 
approach is distributed as the algcurves package with Maple, see [HI [El El]- A purely 
numerical approach to real hyperelliptic Riemann surfaces was given in fT5 s 16j, and for 
general Riemann surfaces in [17J. For a review on computational approaches to Riemann 
surfaces the reader is referred to [3]. 

The codes |11[ HT] compute the periods of a Riemann surface in a homology basis which 
is determined by an algorithm due to Tretkoff and Tretkoff |36j . This homology basis 
is in general not adapted to possible symmetries of the curve (as the involution r of real 
curves). It means that the action on the computed homology basis of any automorphism of 
the curve cannot be expressed in a simple way in terms of this basis. However, the choice of 
a basis, where certain cycles are invariant under the automorphisms, is often convenient in 
applications. In the context of solutions to integrable PDEs on general compact Riemann 
surfaces as for the Kadomtsev-Petviashvili (KP) (see [13J) and the Davey-Stewartson (DS) 
equations (see |3CH I26j). smoothness conditions are formulated conveniently in a homology 
basis adapted to the anti-holomorphic involution r defined on the surface. For instance, 
on a real surface there exists a canonical homology basis (A, B) (that we call for simplicity 
symmetric homology basis in the following) satisfying the conditions (see [35, 38j) 

rA\ = (i g o \ (A 

TB \W -I a \B 



where H £ .M 9 (Z/2Z) is a g x g matrix which depends on the topological type (g,k,a) 
(see section 2); here I g denotes the g x g unit matrix. 
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In |27] we studied for the first time numerically solutions to integrable equations from 
the family of NLS equations, namely, the multi-component nonlinear Schrodinger equation 
and the (2 + l)-dimensional DS equations. A symplectic transformation of the computed 
homology basis to the symmetric homology basis was introduced in |27j and was con- 



structed explicitly for concrete examples. In the present work we give in Theorem 3.3 



a complete description of such a symplectic transformation depending on the topological 
type of the underlying real algebraic curve. The formulas are expressed in terms of period 
matrices of holomorphic differentials on the curve; these holomorphic differentials must 
satisfy the condition T*Uj = i/j, j = l,...,g, where r* is the action of r lifted to the 
space of holomorphic differentials. This allows to give an algorithm to construct explicitly 
a symplectic transformation of an arbitrary canonical homology basis for a real Riemann 
surface to the symmetric form. The algorithm permits to systematically study smooth real 
solutions on general real Riemann surfaces for equations like KP and DS starting from a 
representation of the surface via an algebraic curve which so far was only possible for hy- 
perelliptic surfaces. In addition, this provides a numerical way to compute the topological 
type of a real Riemann surface for given periods of the holomorphic differentials. 

The paper is organized as follows. In section 2 we introduce a symmetric homology basis 
which depends on the topological type of 1Z. In section 3 we give explicitly the symplectic 
transformation between the computed homology basis and the symmetric homology basis. 
This result will be used in section 4 to construct an algorithm which gives the topological 
type (g, k, a) of a real algebraic curve for given periods of the holomorphic differentials. In 
section 5 we discuss examples of real curves for higher genus. Some concluding remarks 
are added in section 6. 



2. Symmetric homology basis 

In what follows 1Z denotes a compact Riemann surface of genus g > 0. A homology 
basis (.4, B) := {A.%, ■ ■ ■ , A g , B\, . . . , B g ) with the following intersection indices 

Ai o Bj = 5{j Ai o Aj = Bi o Bj = 0, 

is called a canonical basis of cycles. With A (resp. B) we denote the vector (Ai, ■ ■ ■ ,A g ) 1 
(resp. (£>i, . . . ,B g ) ). Canonical homology bases are related via a symplectic transforma- 
tion. Let (A, B) and (A, B) be arbitrary canonical homology bases on 1Z. Then there exists 

a symplectic matrix € Sp(2g, Z) such that 



(2.1) 



A B\ (A\ _ (A 
C D [b) ~ \B 



Recall that a symplectic matrix M G Sp(2g,Z) satisfies M t J g M = J g , with the matrix J g 
given by J g = f ^ , where I g denotes the g x g unit matrix. Symplectic matrices 

M = € Sp(2g, Z) are characterized by the following system: 

(2.2) A*D-C t B = I g , 

(2.3) A l C = C l A, 

(2.4) D l B = B l D. 
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Moreover, the inverse matrix M 1 is given by 



(2.5) 



M 



-i 



D* 
-C* 



-B t 
A* 



Now let r be an anti-holomorphic involution defined on 1Z. Recall that (g,k,a) denotes 
the topological type of 1Z, where k is the number of connected components of 1Z(M) (the 
set of fixed points of r), and a = if the curve is dividing (i.e., if TZ \ 1Z(M) has two 
components), a = 1 if it is non-dividing (i.e., if 1Z \ 1Z(M) has just one component). A 
curve with the topological type (g,g + 1,0) is called an M-curve. According to Proposition 
2.2 in Vinnikov's paper [38J, there exists a canonical homology basis (A,B) (called for 
simplicity symmetric homology basis) such that 



(2.6) 



tA 
tB 



where H is a block diagonal g x g matrix which depends on the topological type (g,k,a) 
of 1Z and is defined as follows: 



- if k > and a = 0, 



/0 1 
1 



if k > and a = 1, 



/I 



V 



rank(H) = g + 1 — k in both cases; 

- if k = o, 



/o i 

1 



\ 



1 

1 0/ 



1 

1 



0/ 



0/ 



/o 1 

1 



for even g or 



V 



rank(H) = g if g is even, rank(H) = g — 1 if g is odd. 



\ 



1 

1 



for odd g; 



V 
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Remark 2.1. For given information whether there are real ovals (k > 0) or not (k = 0), 
the matrix H completely encodes the topological type of the real Riemann surface. In [35] 
a different, but equivalent form of HI was used, which has ones only in the antidiagonal, if 
the rank is equal to the genus, or in a parallel to the antidiagonal for smaller rank. 

Example 2.1. Consider the hyperelliptic curve of genus g defined by the equation 

23+2 

(2.7) y 2 = H(x-Xi), 

i=i 

where the branch points x% 6 M are ordered such that x\ < ... < X2 g +2- On such a 
curve, we can define two anti-holomorphic involutions t\ and T2, given respectively by 
Ti(x,y) = (x,y) and T2(x,y) = (x,—y). Projections of real ovals of T\ on the x-plane 
coincide with the intervals [x2 g +2, x\], . . . , [x2g,X2g+l] > whereas projections of real ovals of 
T2 on the x-plane coincide with the intervals [x%, X2], ■ ■ ■ , [x2 g +i, %2g+2\- Hence the curve 



(2.1) is an M-curve with respect to both anti-involutions t\ and T2. 

If all Xi are non-real and pairwise conjugate, the curve has no real ovals with respect to 
the involution T2; it is dividing for the involution t%. 

3. Symplectic transformation between homology bases 

In this section we construct a symplectic transformation between an arbitrary homology 
basis (A,B) and a symmetric homology basis on a real Riemann surface (TZ, r), where r 
denotes the anti-holo morp hic involution defined on 1Z. The main result of this paper is 



contained in Theorem 



3.3 



A B 

which gives the underlying symplectic matrix ( ^ ^ ) in terms 



of the period matrices of holomorphic differentials satisfying the condition (3.2) in the 



homology basis (A,B). The key ingredient in this context is the description of the action 



of r on the cycles (A, B), given in Proposition 3.1 by the integer matrix R. Then Theorem 



states that the column vectors of the matrix (^4 £?) form in fact a "L-basis of the 



integer kernel of the matrix R — I20. The matrix Q in Theorem 3.3 encodes the degree of 



freedom in the choice of such a Z-basis. For the ease of the reader, we start recalling some 



basic facts from the theory of 'L-modules used to prove Theorem 3.3 which differs from 
the usual linear algebra over vector spaces. 

3.1. Basic facts from the theory of Z-modules. A Z-module or more generally a A- 
module where A denotes a commutative ring, is a natural generalization of vector spaces 
where the usual scalar field is replaced by the ring A. For a review on the subject we refer 
to |29j . In what follows we assume that A is the principal ring Z and we denote by A4 a 
Z-module. 

Definition 3.1. 

i. M is of finite type if it admits a finite set of generators. 

ii. M is said to be free if there exists a TL-basis, namely, a set {xi}i^i C A4 with I C N 
such that any element x £ Ai can be written uniquely as x = Yliei a i x i w here the 
scalars ai G Z are non-zero only for a finite number of them. 

The following theorems provide important results in the general theory of modules over 
a principal ring. 



Theorem 3.1. 
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i. If Ai is free and of finite type then all Z-bases of Ai are finite with the same 
cardinality called the rank of Ai. 

ii. Since Z is integral, the rank of Ai equals the dimension of the Q-vector space § -1 .M 
where 8 = Z\ {0}. 

iii. A submodule of a free Z-module of finite type of rank n is a free TL-module of finite 
type of rank r < n. 

The example we will need in this paper are submodules of Z n which are because of the 



Theorem 3.1 free modules of finite type with rank r <n. 



Remark 3.1. Contrary to the case of vector spaces, if Af is a free submodule of the free 
module Ai with the same rank, this does not imply that Af equals Ai. Consider the 
submodule 2Z of Z as an example. 

The following theorem, also called the theorem of the adapted basis, gives the classifica- 
tion of modules over the principal ring Z: 

Theorem 3.2. Let Ai be a free Z-module of rank n and let Af be a submodule of Ai. 
Then there exists a basis (ei, . . . ,e n ) of Ai and unique non-zero integers (pi, . . . ,p r ) (with 
r < n) such that 

(1) (piei, . . . ,p r e r ) is a Z-basis of Af 

(2) pi | j>2 1 • • • \Pr (this means that Vi, pi divides all pj with j > i). 

In what follows GL n (Z) denotes the set of n x n invertible matrices over Z; it is well 
known that the determinant of these matrices equals ±1. Notice that if M denotes the 
matrix formed by the column vectors of a Z-basis of a Z-module of rank n, then other bases 
are given by the column vectors of matrices of the form MQ with Q S GL n (fL). Therefore, 
Theorem |3.2| has the following matrix interpretation that we will use in section |3.3| for 
any non-zero m x n matrix M £ Ai m ,n^f) of rank r, there exist matrices U G GL m (Z), 
V £ GL n (Z) such that 

(3.1) UMV=( D[agiPl Q --- Pr) ° 

where pi E Z\{0} satisfy £>i|p2| • • • \Pr and where Diag(.) denotes the diagonal matrix. This 
is called the Smith normal form of M. In section [4] we will use a well known algorithm to 
compute the Smith normal form. In particular, this algorithm provides a Z-basis of the 



integer kernel of M given by the last n — r column vectors of the matrix V in (3.1 ) 



3.2. Action of r on an arbitrary homology basis. We denote by (A,B) a symmetric 



homology basis (i.e. which satisfies (2.6)). Let (vi, . . . , v g ) be a basis of holomorphic 
differentials such that 



(3.2) T*Vj = uj, j = l,...,g, 

where r* is the action of r lifted to the space of holomorphic differentials: T*uo(p) = u}(rp) 
for any p G 1Z. The matrices P4 and Pq defined by 



(3.3) (P A ) ij= l Uj , (P B ) ij= V j, i,j = l,...,g 

are called the matrices of A- and ^-periods of the differentials Vj. From (2.6) and (3.2) we 
deduce the action of the complex conjugation on the matrices P4 and Pg: 

(3.4) (P A ) tJ G R, 
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(3.5) 



Pb = -Pb + HP4. 



Denote by (A, B) an arbitrary homology basis. From the symplectic transformation 
( |2.1[ ) we obtain the following transformation law between the matrices P A , Pg and P4, Pg 
denned in (3.3): 



(3.6) 



A B 
C D 



P 



A 



Pi 



Pa 
Pb 



Therefore, by (3.4) one gets 

(3.7) 
(3.8) 



ARe(P A ) +BRe(P B ) = P A 
Alm(P A ) +Blm(P B ) = 0, 



and by (3.5 ) 

(3.9) 
(3.10) 



CRe(P A )+DRe(P B ) = -MP A 
C Im (P A ) + D Im (Pg) = Im (P B ) . 



From (3.7) and (3.10) it can be checked that the matrices ^4Re {P A ) + PRe (Pg) and 
C Im (P4J + D Im (Pgj are invertible (the first because P4 is, for the second see [27] for 
more details). The following Lemma proved in |27j shows that it is sufficient to know the 
pairs of matrices A, B or C, D to get the full symplectic transformation (3.6): 



Lemma 3.1. The matrices A,B,C,D G M g (Z) solving (3.7)-(3.10) satisfy: 

(3.11) A t = Im(P B ) [CIm(P A )+DIm(P B )]' 1 

(3.12) B t = -Im{P A ) [C Im (P A ) + D Im (Pg)] 1 

(3.13) C* 

(3.14) D t 



^A f M- Re(P 8 ) [ARe(P A ) +BRe(P 8 )] 1 
\ P*H + Re (P A ) [A Re (P A ) + B Re (Pg)] 1 



The action of r on an arbitrary homology basis (A, B) can be written as 



(3.15) 



tA 
tB 



R 



where R G M2g{^)- I n the following proposition, we give an explicit expression for the 
matrix R in terms of the period matrices P A and Pg only. 



Proposition 3.1. The matrix R defined in (3.15) is given by 
(3.16) 



R 



2Re(P, 



B) 



1 Im(p t 7 ) +1 



2Re(P, 



B) 



A 

- 1 Im ( Pi 



-2Re(P 



A) 



2Re(P, 



B) 



A- 1 Im (Pjj 
1 Im (P A ) + I t 



where 
(3.17) 



Im (Pi) Re (P A ) - Im (pU Re (Pg) . 
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Proof. Using (2.6) we deduce the action of r on (2.1): 

-i 



(3.18) 

which yields 
(3.19) 



A B 
C D 



R 



I, 



A B 
C D 



R 



A B 
C D 



A B 
C D 



In other words, a symplectic matrix which transforms a basis (A, B) to a symmetric form 



also satisfies (|3.19|). From (|3.19|) and (|2.5|) one gets 
(3.20) 



R 



2 C l B - A l WB + l g y 
-2 C l A + A*MA 



2D l B - B l UB 
(2C l B - A l WB + I g )J ' 



Replacing A via (3.11) and B via (p. 12) in (3.13) and then using (3.11) again to eliminate 



the factor [C Im ( P4) + D Im (Pg)J leads to a relation for C l A only (using in the last step 
(3.12) instead of (3.11) gives a relation for C l B). Similarly one gets a relation for D l B, 



2 C l A = A l WA - 2 Re (P B ) M" 1 Im ( P { 
2 C f B = A l UB + 2 Re (Pg) MT 1 Im [p 1 ^ 
2 D f B = B f MB - 2 Re (P A ) MT 1 Im (p*A , 



with M given by (3.17). Substituting these relations in (3.20) one gets (3.16). 



□ 



3.3. Symplectic transformation. In this part we present in Theorem 3.3 the main result 
of the present paper: the symplectic transformation between an arbitrary basis (A, B) on 
a real Riemann surface and a homology basis adapted to the symmetry is given in terms 
of the period matrices Pi,Pg defined in the previous section. This result will allow to 
construct in section [4] an algorithm which computes the topological type of a given real 
Riemann surface. 

We start with the following lemma which describes the spectral properties of the matrix 



R (3.16): 



Lemma 3.2. The matrixH in (3.16) is diagonalizable over Q with eigenvalues 1 and — 1. 



The dimension of the corresponding eigenspaces equals g. 
Proof. It is straightforward to see that the matrix 



is diagonalizable over 



the eigenvalues are 1 and —1, and the dimension of the corresponding eigenspaces equals 
g. Therefore, by (3.19) the same holds for the matrix R. □ 



In what follows we denote by 
(3.21) lC z :={w£Z 2 V;(R t -h g )w = 0} 

the integer kernel of the matrix R* — l2 g . According to the theory of modules over the 
principal ring Z, the Z-module K,% admits a Z-basis which can, for instance, be computed 



from the Smith normal form (see section 3.1). 
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Theorem 3.3. Let TZ be a real compact Riemann surface of genus g and (A, B) a canonical 
homology basis on TZ. For the given matrices of periods P%,Pg of holomorphic differentials 



satisfying (3.2), let 



be the 2g x g matrix formed by a Z-basis of K.% (3.21). Then a 



symplectic matrix in (2.1) which relates the homology basis (A, B) to a symmetric homology 

Q 



basis (A, B) is given by: 
' At \ _ {Si 



(3.22) 
(3.23) 



R<(P D H)[S[Re(P A ) + SiRe(P,)] 



Re(P 



A 



-1 



(Q) 



t\-i 



where the matrix H is defined in section 2, and where Q G GL g (Z) is such that 



(3.24) 



QMQ 1 



-Re (Pp 
Re(P 



A, 



P{ [Si Re (P A ) + S f 2 Re (P $ )] (mod 2). 



Proof. Notice that (3.19) can be rewritten as 



which, in particular, gives the following condition for the matrices A and B: 
(3.26) (A B)(R-h g ) = 0. 



(A B\ 

Denote by u\, . . . , u g , i>i, . . . , v g the column vectors of the matrix ( ^ Jj) ' ^ (3.26) 



the vectors U{ for i = 1, ... ,5 lie in the integer kernel Kz of the matrix R* — I20. Let us 



prove that (u\, . . . ,u g ) form in fact a Z-basis of the module ICz- By Lemma 3.2 one has 
dim(K,Q) = g, where JCq denotes the kernel of R* — l2 g over the field Q, which by Theorem 



3.1 yields rank()Cz) = g. Therefore one has to check that (u\, . . . , u g ) are free vectors over 



Z and generate the Z-module /C^. Notice that here it is important to check that these 



vectors form a set of generators for the Z-module /C^, as we saw in Remark 3.1 
(A B\ 1 

Since ( ^ ^1 G GL2 9 (Z), the 2g vectors u\, . . . ,u g ,v\, . . . ,v g form a Z-basis of the 

module Z 2ff , which in particular implies that these vectors are free over Z. Then it remains 
to prove that the vectors Ui,i = l,...,g generate the Z-module fCz which is done by 
contradiction as follows. Let w G /C^ which we write in Z 2s as w = 5^f=i a i u i + Sj=i Pj v j 
with cti, f3j G Z such that at least one of the f3j is non-zero. Since w, ui, . . . , u g G /Cz, one 
has v := Ylj=i Pj v j ^ ^-z an( i v 7^ 0- We deduce that m, . . . , u g ,v are g + 1 free vectors 
in ICz- This is impossible since rank(ICz) = g- Thus one has f3j = for j = 1, . . . , g which 
implies that the vectors Ui,i = 1, . . . , g generate the Z-module KL%. 
Hence we can write 

(3.27) (|) = (|) Q, 

for some Q G GL 9 (Z), where the g column vectors of the matrix form a Z-basis 

of the module tC-£- Here the matrix Q encodes the freedom in the choice of such a basis. 



The matrices C and D are then given by (3.13) and (3.14). It follows that these two 
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matrices are integer matrices if and only if the matrix Q satisfies (3.24), which completes 
the proof. □ 



Remark 3.2. If the curve is an M-curve, then 
Theorem 3.3 are given by 



and the matrices A, B, C, D in 



(3.28) 
(3.29) 



A* 
B f 



Q 



-Re (Pg) 



- 1 fnfi\-l 



[^ 1 Re(P | )+^ 2 Re(Pg)] (Q 



Re(P A ), 
where Q G GL g (X) is arbitrary. 

Remark 3.3. As explained below, if the curve is not an M-curve, namely, HI 7^ 0, one can 
construct explicitly a matrix Q G GL g (Z/2Z) such that (3.24) holds. 



This construction of Q is based on the Smith normal form of the matrix 



section 



sec 



3.1) which allows the simplification of the system (3.24). 



Lemma 3.3. There exist U G GL^g^E), V G G7 g (Z) such that 

'Si 



(3.30) 



U 



V 



where £ is a g x g diagonal matrix with elements £ a = ±1 for i = 1, . 



,9- 



Proof. By (|3TT]) there exist U G GL 2g (Z), V G GL g (Z) and pi, . . . ,p g G N \ {0} satisfying 
Pi |P2 1 - - - \Pg such that 



(3.31) 



U 



V 



where D := Diag(pi, . . . ,p g ). The fact that 
equalities: 

'Si C l 



£ can be deduced from the following 



(3.32) 1 = det p 



det 



S 2 D* 



det(Q) = ±det 



F x 
F 2 



±det(D)det(F 2 ) 



where we multiplied the matrix in the determinant from the left by U and from the right by 
^ ^ , and used det(Q) det (17) det(V) = ±1 since Q,V G G7 9 (Z) and [/ G GL 2g (Z); 

^ j = [/ ( jj t J . We deduce that det(D) = ±1 since the right-hand side of (3.32 ) is a 
product of determinants of matrices with integer coefficients. This completes the proof. □ 

Now let us define matrices Ni,N 2 G A4 g (Z) as follows: 



here 



(3.33) 

Then one has: 



Ni 
N 2 



2VU 



-Re [P f . 
Re(P 



A 



?_ {)[SiRe(P A )+SlRe(P s )]-\ 



Proposition 3.2. Q G GLgiT) satisfies (3.24) if an d only if it solves 
(3.34) QHQ* = Ni (mod 2), 
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where N% G M g (Z) is defined in (3.33). Moreover, if Q is a particular solution of (3.34) 
then the general solution can be written as QQq where Qq solves 



(3.35) 



QoHQo = H (mod 2) 



Proof. Multiplying equality (3.24) from the left by the matrix U of Lemma 3.3 and using 



(3.30 ) one gets 



V~ 1 QWQ t = 2U 



-RefP, 



Re(P 



A 



which is equivalent to 
'I 



QM.Q 1 = 2VU 



-Re (P B 



Re IP 



A 



[SlReiPjd+SlReiPg)] 1 (mod 2), 



[S[ Re (P|) + S\ Re (Pg) ] 1 (mod 2) . 



Using the definition (3.33) one gets d3.34). 



Now to check that a transformation of the form Q — > QQq where Qq solves (3.35) is 



the only one which preserves (3.34), notice that such a transformation corresponds to a 



symplectic transformation between the symmetric homology basis obtained from Theorem 



3.3 and another symmetric homology basis, since this coincides with a change of the Z- 



basis in (3.22). Hence from (2.6) it is straightforward to see that the symplectic matrix 



which relates two symmetric homology bases is given by 
(3.36) 



Qo 





1-1 



with Qq £ GL g (Z) satisfying QqMQ 



KWQl-Q^W) Q 

M (mod 2). This completes the proof. 



□ 



4. Algorithm for the computation of the topological type (g, k, a) 



The results of the previous section allow us to formulate an algorithm to transform an 
arbitrary canonical homology basis (A,B), for instance obtained via the algorithm [36], to 
a symmetric basis {A,B) satisfying (2.6). The key task in this context is the computation 
of the matrix Q in ( 3.34[ ). In the process of computing Q, the matrix EI giving the topology 
of the real Riemann surface can be determined. 

The starting point of the algorithm are the periods P% and Pg of a basis of differentials 



u\, . . . , v g ) satisfying (3.2 ). Notice that condition (3.2 ) is important. The Maple algcurves 



package generates such differentials for real curves by default. For the Matlab code used 
in this paper this is in general not the case if a numerically optimal approach is used 
to determine the holomorphic differentials, see [17J for details. However, it is possible 
to determine a basis of the holomorphic differentials with rational coefficients which will 



satisfy condition (3.2). For the examples in the following section, we always choose this 
option. 



With these periods the code computes the matrix R via (3.16). This matrix will have 
integer entries up to the used precision (by default 10 -6 in Maple and 10 -12 in Matlab). 
Rounding has to be used to obtain an integer matrix. Computing the Smith normal form 



(3.1) of the matrix R* — J.2gi we obtain from the last g vectors of the resulting matrix V 

Si 

S-2 



(3.27). 



a Z-basis of the integer kernel K% (3.21), i.e., the column vectors of the matrix 



m 



An algorithm to compute the Smith normal form UMV = S of an integer matrix M 
is implemented in Maple. This algorithm can be called from Matlab via the symbolic 
toolbox. We use here an own implementation of the standard algorithm to compute the 
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Smith normal form which we briefly summarize: we always work on column j starting 
with j = 1. If there is no nonzero element in this column, it is swapped by multiplication 
with an appropriate matrix U from the left with the last column with a nonzero element. 
If the element Mjj = 0, a row with a nonzero element in position j is added (to avoid 
clumsy notation, the transformed matrix is still called M). Then all nonzero elements Mj^ 
for k ^ j are eliminated by adding row j with appropriate multipliers obtained via the 
Euclidean algorithm. In the same way the row with index j is cleaned by acting on M via 
multiplication by a matrix V from the right. If the resulting element Mjj does not divide 
all other elements of M, one of these elements is added by using the Euclidean algorithm 
to Mjj in a way that the latter becomes smaller. This destroys possibly the nullity of 
the remaining elements in column j and row j which thus have to be cleaned as before. 
This process is repeated until Mjj divides all other elements of M. Then the index j is 
incremented by 1. The procedure is repeated until the Smith normal form is obtained. 

Notice that this standard algorithm has a well known problem: in general for larger 
matrices the entries of the matrices U and V become very large. Though these are integer 
matrices, this is problematic once some of the entries are of the order of 10 16 (machine 
precision in Matlab is 10~ 16 which implies that integers of the order of 10 16 , which are 
internally treated as floating point numbers, can no longer be numerically distinguished). 
There are more sophisticated algorithms to treat larger matrices as the one given in [22J. 
In practice the standard algorithm works well for examples of a genus g < 6 which is 
sufficient for our purposes. Only for higher genus, the algorithm [22 j w ould be needed. 

By computing the Smith normal form of the vector (^gj ^ n (3-30), we get the needed 

quantities to compute the matrix Aq in (3.33). The main task is then to determine the 



equations (3.22 



matrix Q i n ( 3.2 4 ) sin ce for given periods, the whole symplectic matrix (2.1) follows from 

'Si 



and (3.23) for given 



S, 



The matrices Q and 



from relation (3.34) for a given matrix Aq by standard Gaussian elimination in 
by imposing the block diagonal form of section 2 on EI as we will outline below. 



can be determined 
/2Z and 



Remark 4.1. It was shown in section 2 that the matrix 
diagonal or to consist of blocks of the form 



can be chosen to be either 



1 

1 



A unique determination of the matrix H in the computation of Q from (3.34) is only 



possible, if this block Hq cannot be related through a similarity transformation in Z/2Z to 
I2. In fact this is the case since Hq cannot be diagonalized in Z/2Z. The same reasoning 
applies if the matrix H consists of several blocks Hq and zeros otherwise. 
However, a block of the form 

/l 0^ 
Hi = 1 

V° 1 °y 

can be diagonalized in Z/2Z by multiplication from the left and the right by a matrix of 
the form 

(I 1 1 N 
Qi = 1 1 
U l 0, 
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It follows that if the matrix H determined in the computation of Q from (3.34) has a 
non-zero diagonal element, then 



has the form H below (H, 



can be diagonalized in several steps: if the matrix 

k,k + l) 



/I 



H 



5ij for i,j < k and a block Hq for i, j 
\ /111 



1 

1 



Q 



\ 



1 

1 



V / V / 

then multiplication from the left and the right with a matrix Q (not shown elements of this 
matrix are 0) gives the identity matrix for i,j < k + 1. Applying this procedure several 
times will lead to a diagonal matrix HI. 



The algorithm for the computation of Q and H via the similarity relation ( |3.34 ) for given 



Ni by imposing a block diagonal form for the matrix HI as in section 2 uses in principle 
standard Gauss elimination on the rows and columns of N\ over the field Z/2Z with minor 
modifications as detailed below. We only describe the action on the columns via a matrix 
q from the right, since the action on the rows follows by symmetry by multiplication with 
q t from the left: 

- if Ni = (mod 2), put HI = and Q = I g and end the algorithm, otherwise, put the index 
j of the column under consideration equal to 1; 

- if column j contains only zeros, it is swapped with the last column with non-zero entries; 

- if there is a 1 in position j of the column, all further non-zero entries in the column are 
eliminated in standard way; 

- if there is a 1 in the column, but not in position j, rows are swapped in a way that it 
appears in the position j + 1 of the column (it cannot be put to position j as explained in 



Remark 4.1). Further ones in the column are eliminated; 



- if there was a non-zero diagonal entry, the column index j is incremented by 1, if there 
was a block ^ ^ , the index j is incremented by 2. Then the algorithm is repeated with 
column j until j = g or until the columns wi th in dex j and higher only contain zeros; 



if there are blocks of the form H in Remark 4.1 then HI wi ll be diagonalized by multipli- 



cation with the corresponding matrix Q as explained in 4.1 



5. Examples 

In this section we study examples of real algebraic curves, provide the computed periods 
and the application of the algorithm to obtain a symmetric homology basis as well as the 
matrix HI encoding the topological information of the curve. For convenience we use here 
the Matlab algebraic curves package, but the same examples can be of course studied with 
the Maple package. We also give graphical representations of the real variety of an algebraic 
curve if there is any which is generated via contour plots of f(x, y) = for real x and y (this 
corresponds to the command plot_real_curve in the Maple algcurves package). This is 
not identical with the set 7£(R) of real ovals of the Riemann surface since the curves may 
have singularities, whereas the Riemann surface is defined by desingularized such curves 
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(see [17J for how this is done in the Matlab package). Thus there may be cusps and self- 
intersections in the shown plots. Moreover these plots are not conclusive if the curves 
come very close, and if there are self-intersections as in Fig. [4] below. In addition we only 
show the curves for finite values of x and y from which it cannot be decided which curves 
cross at infinity and which lines belong to the same ovals as in Fig. [5] They just serve 
for illustration purpose, for more sophisticated approaches, see [2j |33l HU H21 123 134] , 
The computed number of real ovals via the algorithm is, however, unique: as outlined in 
section 2, it follows from the rank of the matrix H (for k ^ one has k = g + 1 — rank(H)) . 
We always assume in the following that it is known whether there are any real ovals. This 
allows the unique identification of the topological type (g,k,a) via the matrix HI. 
The Trott curve |37] given by the algebraic equation 

(5.1) 144 (x 4 + y 4 ) - 225 (x 2 + y 2 ) + 350 x 2 y 2 + 81 = 

is known to be an M-curve of genus 3 (it has the maximal number g + 1 = 4 of real ovals, 
as can be seen in Fig. [TJ. Moreover, this curve has real branch points only (and 28 real 
bitangents, namely, tangent lines to the curve in two places). Our computed matrices of 



FIGURE 1. Real ovals of the Trott curve (5.1), an M-curve of genus 3. 



A and i3-periods denoted by aper and bper respectively reacj^] 
aper = 

-0.0000 + 0.0235i -0.0000 + 0.0138i -0.0000 + 0.0138i 
+ O.OOOOi 0.0000 + 0.0277i + O.OOOOi 

-0.0315 -0.0000 + O.OOOOi 0.0250 - O.OOOOi 



bper 



^For the ease of representation we only give 4 digits here, though at least 12 digits are known for these 
quantities. 
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-0.0315 + 0.0235i 
-0.0000 + O.OOOOi 
-0.0000 - 0.0235i 



-0.0000 + 0.0138i 
-0.0250 + 0.0277i 
0.0000 + 0.0138i 



-0.0250 + 0.0138i 
0.0250 - O.OOOOi 
+ 0.0138i. 



For this the algorithm produces as expected 



and Q of (3.34) the identity matrix. 



The symplectic transformation found by the algorithm via (3.22) and (3.23) has the form 
[A.B.C.D] = 



1 








-1 





-1 












1 















1. 



Since we will not actually use the matrices A, B, C, D in this article, and since they follow 
for given periods from H and Q via (3.22) and (3.23) we will only give them for this 
example. However they are needed e.g. for the study of algebro-geometric solutions to 
integrable equations as KP and DS as in [27J. 
The Klein curve given by the equation 



(5.2) 



y 



x(x — 1) 







has the maximal number of automorphisms (168) of a genus 3 curve. The computed periods 
read 



aper 



-0.9667 + 0.7709i 
-1.2054 - 0.2751i 
-0.4302 - 0.8933i 



0.9667 + 0.2206i 
-0.4302 + 0.8933i 
1.7419 + 1.3891i 



0.9667 - 2.0073i 
-1.7419 + 1.3891i 
-1.2054 + 0.2751i 



bper = 

-2.7085 - 0.6182i 
-2.1721 - 1.7322i 
0.9667 + 0.2206i 

The algorithm finds 



-0.2387 + 0.4958i 
0.5365 - 0.1224i 
-0.9667 + 2.0073i 



1.3969 - 1.1140i 
-0.7752 - 1.6097i 
-0.9667 + 0.7709i. 



1 











1 











1 



111 

1 
1 0. 

Therefore, the topological type of the curve is (3, 1, 1), namely, the curve has genus 3, one 
real oval (as can be also seen in Fig. [2]) and is non-dividing. 
The Fermat curve 

(5.3) y n + x n + 1 = 
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FIGURE 2. Real variety of the Klein curve (5.2 ), the curve with the maximal 
number of automorphisms in genus 3. 



has for n = 4 the topological type (3,0, 1). This is confirmed by the algorithm. For 
periods 

aper = 






9270 


+ 





OOOOi 





0000 







9270i 





0000 


- 


9270i 


-0 


0000 


+ 





OOOOi 





0000 


+ 





OOOOi 





0000 


- 1 


854 li 







+ 





9270i 


-0 


9270 


+ 





OOOOi 





0000 


- 


9270i 


bper 































9270 


+ 





9270i 





9270 







9270i 





0000 


+ 


OOOOi 





0000 








-0 


9270 


+ 





9270i 





9270 


- 


9270i 


-0 


9270 


+ 





OOOOi 





0000 







9270i 





0000 


- 


9270i 



we find 
H = 






1 





1 


















10 
1 
10 

in accordance with the expectation. For n = 5 the curve has genus 6. We find with 
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aper = 

Columns 1 through 4 

-0.1623 + 0.4995i 
-0.3246 - 0.00001 
-0.4249 + 0.3087i 
-0.4249 - 0.3087i 
0.1003 - 0.30871 
0.5252 + 0.00001 

Columns 5 through 6 

0.4948 + 0.3595i 
0.9896 - 0.00001 
0.1890 - 0.5817i 
0.1890 + 0.58171 
0.8006 + 0.5817i 
0.6116 - O.OOOOi 

bper = 

Columns 1 through 4 

-0.6875 + 0.4995i 
0.1003 - 0.3087i 

-0.6875 - 0.49951 

-0.1623 - 0.4995i 
0.4249 - 0.3087i 

-0.1623 + 0.4995i 

Columns 5 through 6 

-0.1168 + 0.3595i 
0.8006 + 0.58171 
-0.1168 - 0.3595i 
0.4948 - 0.3595i 
-0.1890 + 0.58171 
0.4948 + 0.3595i 
the matrices 
H = 



1 











1 











1 






























0. 


1890 


+ 


0. 


,5817i 


0. 


3780 


+ 


0. 


,00001 


0. 


0722 




0. 


2222i 


0. 


1890 




0. 


, 1373i 


0. 


3058 


+ 


0. 


,2222i 


u . 






u . 


, ooyoi 


0. 


4249 


- 


0. 


,3087i 


0. 


8498 


- 


0. 


OOOOi 


1. 


1125 


+ 


0. 


,8082i 


0. 


4249 


- 


1. 


,30781 


0. 


2626 




0. 


,8082i 


0. 


1623 


- 


0. 


,4995i 


0. 


8006 


+ 


0. 


,5817i 


0. 


1168 


- 


0. 


,3595i 


0. 


3058 




0. 


2222i 


0. 


2336 


+ 


0. 


,00001 


0. 


6116 




0. 


,0000i 


0. 


4948 


- 


0. 


,3595i 


0. 


1003 




. 


OA07 • 

, 3087i 


0. 


6875 




. 


, 4995i 


0. 


2626 


+ 


. 


, 8082i 


1. 


3751 




0. 


OOOOi 


0. 


5252 


+ 


0. 


OOOOi 


0. 


1623 




0. 


,49951 


c 


i 
























i 










1 















1 















1 



-0. 


, 1623 


+ 





.4995i 


-0. 


,3246 


+ 





.00001 


0. 


,5252 







.OOOOi 


0. 


, 1623 


+ 





,4995i 


0. 


,3246 








0. 


,1623 


+ 





,4995i 


-0. 


,6875 


+ 





.4995i 


0. 


,1003 







.3087i 


0. 


,2626 


+ 





.80821 


-0. 


, 1623 


+ 





,4995i 


0. 


,4249 


+ 





.30871 


-0. 


,5252 












.4948 


+ 


0. 


,3595i 





.9896 




0. 


OOOOi 


0. 


.6116 




0. 


OOOOi 


-0 


.4948 


+ 


0. 


,3595i 


-0 


.9896 




0. 


OOOOi 


-0 


.4948 


+ 


0. 


,3595i 


-0 


.1168 


+ 


0. 


,3595i 





.8006 


+ 


0. 


,5817i 


0. 


.3058 




0. 


,2222i 


0. 


.4948 


+ 


0. 


,3595i 


-0 


.1890 




0. 


,5817i 


-0 


,6116 
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q = 



1 

















1 








1 




















1 











1 





1 








1 


1 





1 











1 





1 


1 



This implies that there is one real oval as can be also seen in Fig. [3] and the curve is 
non-dividing. This corresponds to (g, k, a) = (6, 1, 1). 




FIGURE 3. Real variety of the Fermat curve (5.3) for n = 5 



For the curve 
(5.4) 

of genus 3 we have 
aper = 



y 3 - 2x 3 y - x 9 = 






4021 - 





6964i 


-0 


6748 - 


- 1 


1688i 





5985 


+ 


1 


0367i 





4764 - 





9006i 


1 


5026 - 


- 


1823i 


2 


0418 


+ 





5711i 


1 


5598 - 





9006i 





3157 - 


- 


1823i 


-0 


9892 


+ 





5711i 



bper 



1 


1577 


+ 


2041i 





3591 


- 


9865i 





3907 


+ 


4656i 





5417 


- 


3128i 





5934 


+ 


3426i 


1 


5155 


+ 


8750i 





6160 


+ 


1086i 


-0 


2344 


+ 


6439i 


-1 


1249 


- 1 


3406i 



which leads to 
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H = 



1 











1 











1 



q = 

ill 

1 
1 0. 

This gives the topological type (3, 1, 1). In particular, the number of real ovals equals one. 
The real variety of the curve, which has a self intersection and a cusp, can be seen in Fig. H] 




-1.5 -1 



FIGURE 4. Real variety of the curve (5.4) 



The curve 
(5.5) 

-I80x 5 +396yx 4 -307x 3 y 2 +107x 2 y 3 +273x 3 -318x 2 y-17xy i +ll7xy 2 -68x+y 5 -l2y 3 +19y 

is known (see |12| ) to be a dividing curve of genus 6. In fact we get for 
aper = 

Columns 1 through 4 



0.0414 + 0.0278i 

-0.1149 - 0.0446i 

0.1149 + 0.0169i 

0.0000 - O.Ollli 

0.0820 - 0.0278i 



-0.0345 + 0.0272i 

0.0000 - O.OOOOi 

-0.0000 + 0.0544i 

0.0000 + 0.0183i 

0.0000 + O.OOOOi 



-0.0979 + 0.0264i 
0.0532 - 0.0226i 

-0.0532 - 0.0805i 
0.0000 - 0.0303i 

-0.0527 - 0.1031i 



-0.3041 + 0.0603i 

0.0000 - O.OOOOi 

-0.0000 + 0.1206i 

0.0000 + 0.1114i 

0.0000 + O.OOOOi 
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-0.0820 - O.OOOOi -0.0000 - 0.0544i 0.0527 - O.OOOOi 



Columns 5 through 6 



-0. 


.5369 


+ 


0. 


0872i 


-2. 


,8149 


+ 0. 


,5433i 





,2427 


+ 


0. 


.00991 


0. 


.7641 


+ 0. 


,03571 


-0. 


.2427 




0. 


,2372i 


-0. 


,7641 


- 0. 


,5381i 





,0000 




0. 


. 1843i 


0. 


,0000 


- 1. 


, 1224i 


-0. 


.0881 




0. 


,2274i 


-0. 


,1275 


- 0. 


,5024i 





.0881 




0. 


,00001 


0. 


,1275 


- 0. 


,00001 



bper = 

Columns 1 through 4 






,0414 


+ 


0. 


,0278i 


-0. 


.0345 


- 0, 


,0094i 


-0. 


0979 


+ 0, 


,0264i 


-0. 


.0089 




0. 


.10091 


-0. 


.0666 


+ 0. 


,01801 


0. 


1091 


- 0. 


.07331 


-0 


,0089 




0. 


,0563i 


-0. 


.0666 


- 0. 


,0180i 


0. 


1091 


- 


.0507i 


0, 


.0320 




0. 


,0000i 


0. 


.0000 


- 0. 


0183i 


0. 


1425 


- 0. 


.0000i 





.1060 


+ 


0. 


.02861 


0. 


,0666 


+ 0, 


,0724i 


0. 


0559 


- 


,0525i 


-0. 


.0580 


+ 


0. 


.01601 


0. 


,0666 


- 0, 


,03631 


0. 


1614 


+ 0. 


.07511 



Columns 5 through 6 



0, 


.5369 


+ 


0. 


0872i 


-2. 


,8149 


+ 


0. 


,5433i 


0, 


.3380 




0. 


1055i 


0. 


,9555 




0. 


,2619i 


0, 


.3380 




0. 


1154i 


0. 


,9555 




0. 


,29761 


0. 


.8311 




0. 


00001 


4. 


,8657 




0. 


,00001 


0, 


.0954 




0. 


1120i 


0. 


, 1914 




0. 


,2048i 





,2716 


+ 


0. 


1021i 


0. 


,4464 


+ 


0. 


, 1691i 



the matrices 



-0, 


,0000 


- 0. 


1206i 


-0 


.3041 


- 0. 


1626i 


-0 


.1162 


+ 0. 


00461 


-0 


.1162 


- 0. 


0046i 





.0000 


- 0. 


1114i 


0. 


.1162 


+ 0. 


1252i 





.1162 


- 0. 


11601 



H 




1 







1 













1 









1 

























Q 



1 



1 








1 
1 
1 








1 











1 







1 
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This means that it is a dividing curve with 3 real ovals, which is, however, not obvious 
from Fig. [5] The topological type is (6, 3, 0). 




FIGURE 5. Real variety of the curve (5.5) 



6. Outlook 

The algorithm presented in this paper allows the transformation of an arbitrary canonical 
homology basis to a form adapted to the underlying symmetry, here the anti-holomorphic 



involution. This permits to find a basis satisfying relation (2.6). A similar condition can 



be imposed for any involution, and the algorithm presented here can be easily adapted to 
that case. 

In general, the presence of symmetries allows to significantly simplify the Riemann 
matrix of a surface, but only in a homology basis adapted to the symmetries, for instance 
in a basis such that the ^-cycles are invariant under symmetry operations, see [3] and 
[5] for the Klein curve. The latter reference uses an approach to this problem based on 
Comessati's theorem [7] via two pieces of software, extcurves and CyclePaintei^ It will 
be the subject of further work to generalize the algorithm to symmetry groups beyond 
involutions. In a first step it would be interesting to extend the approach presented in this 
paper to automorphisms r satisfying r ra = id with n > 2. 



References 

[1] V. I. Arnold, The situation of ovals of real plane algebraic curves, the involutions of four- dimensional 
smooth manifolds, and the arithmetic of integral quadratic forms, Funkcional. Anal, i Prilozen 5, 1-9 
(1971). 

[2] D. S. Arnon and S. McCallum, A polynomial time algorithm for the topological type of a real algebraic 

curve, Journal of Symbolic Computation, 5:213-236 (1988). 
[3] E. Belokolos, A. Bobenko, V. Enolskii, A. Its, V. Matveev, Algebro-geometric approach to nonlinear 

integrable equations, Springer Series in nonlinear dynamics (1994). 



2 Located at http://gitorious.org/riemanncycles 



22 



C. KALLA AND C. KLEIN 



[4] A.I. Bobenko, C. Klein, (ed.), Computational Approach to Riemann Surfaces, Lect. Notes Math. 2013 
(2011). 

[5] H. W. Braden, T.P. Northover, Klein's Curve, J. Phys. A 43, 434009 (2010). 

[6] H. Braden, V. Enolskii, T. Northover, Maple packages extcurves and CyclePainter, available at gito- 
rious.org. 

[7] A. Comessati, Sulla connessione delle superficie algebriche reali, Annali di Mat. (3) 23, 215-283 (1915). 
[8] M. Coste, M.-F. Roy, Thorn's lemma, the coding of real algebraic numbers and the computation of the 

topology of semi-algebraic sets, Journal of Symbolic Computation, 5:121-129 (1988). 
[9] B. Deconinck, M. van Hoeij, Computing Riemann matrices of algebraic curves, Physica D, 28, 152-153 

(2001). 

[10] B. Deconinck, M. Heil, A. Bobenko, M. van Hoeij, M. Schmies, Computing Riemann Theta Functions, 
Mathematics of Computation 73, 1417 (2004). 

[11] B. Deconinck, M. Patterson, Computing with plane algebraic curves and Riemann surfaces: the algo- 
rithms of the Maple package "algcurves", in A.I. Bobenko, C. Klein, (ed.), Computational Approach 
to Riemann Surfaces, Lect. Notes Math. 2013 (2011). 

[12] B.A. Dubrovin, Matrix finite-zone operators, Revs. Sci. Tech. 23, 20-50 (1983). 

[13] B. Dubrovin, S. Natanzon, Real theta-function solutions of the Kadomtsev-Petviashvili equation, Math. 

USSR Ivestiya, 32:2, 269-288 (1989). 
[14] H. Feng, Decomposition and Computation of the Topology of Plane Real Algebraic Curves, Ph.D. 

thesis, The Royal Institute of Technology, Stockholm (1992). 
[15] J. Frauendiener, C. Klein, Hyperelliptic theta functions and spectral methods, J. Comp. Appl. Math. 

(2004). 

[16] J. Frauendiener, C. Klein, Hyperelliptic theta functions and spectral methods: KdV and KP solutions, 

Lett. Math. Phys., Vol. 76, 249-267 (2006). 
[17] J. Frauendiener, C. Klein, Algebraic curves and Riemann surfaces in Matlab, in A. Bobenko and 

C. Klein (ed.), Riemann Surfaces-Computational Approaches, Lecture Notes in Mathematics Vol. 

2013 (Springer) (2011). 

[18] A. Gabard, Sur la topologie et la geometrie des courbes algebriques reelles, Ph.D. Thesis (2004). 

[19] L. Gonzalez- Vega, I. Necula, Efficient topology determination of implicitly defined algebraic plane 

curves, Computer Aided Geometric Design, 19:719-743 (2002). 
[20] D.A. Gudkov, Complete topological classification of the disposition of ovals of a sixth order curve in 

the projective plane, Gor'kov. Gos. Univ. Ucen. Zap. Vyp., 87:118-153 (1969). 
[21] A. Harnack, Ueber die Vieltheiligkeit der ebenen algebraischen Curven, Math. Ann., 10, 189-199 

(1876). 

[22] J. Hafner and K. McCurley, Asymptotically fast triangularization of matrices over rings, SIAM J. 

Comput. 20, 1068-1083 (1991). 
[23] D. Hilbert, Uber die reellen Ziige algebraischen Kurven, Math. Ann., 38, 115-138 (1891). 
[24] D. Hilbert, Mathematische Probleme, Arch. Math. Phys., 1:43-63, (German) (1901). 
[25] H. Hong, An efficient method for analyzing the topology of plane real algebraic curves, Mathematics 

and Computers in Simulation, 42:571-582 (1996). 
[26] C. Kalla, New degeneration of Fay's identity and its application to mtegrable systems, preprint 

|arXiv:1104.2568y i (2011). 
[27] C. Kalla, C. Klein, On the numerical evaluation of algebro-geometric solutions to integrable equations, 

Nonlinearity 25 569-596 (2012). 
[28] F. Klein, On Riemann's theory of algebraic functions and their integrals, Dover (1963). 
[29] S. Lang, Algebra, Revised Third Edition, Springer- Verlag, New York (2002). 

[30] T. Malanyuk, Finite-gap solutions of the Davey-Stewartson equations, J. Nonlinear Sci, 4, No. 1, 1-21 
(1994). 

[31] I. Petrovsky, On the topology of real plane algebraic curves, Ann. Math., 39, No. 1, 187-209 (1938). 

[32] K. Rohn, Die Maximalzahl und Anordnung der Ovale bei der ebenen Kurve 6. Ordnung und bei der 
Flache 4. Ordnung, Math. Ann., 73, 177-229 (1913). 

[33] T. Sakkalis, The topological configuration of a real algebraic curve, Bulletin of the Australian Mathe- 
matical Society, 43:37-50 (1991). 

[34] R. Seidel, N. Wolpert, On the exact computation of the topology of real algebraic curves, In Proc 21st 
ACM Symposium on Computational Geometry, p 107-115 (2005). 



23 



[35] M. Seppala and R. Silhol, Moduli Spaces for Real Algebraic Curves and Real Abelian Varieties, Math. 
Z. 201, 151-165 (1989). 

[36] C.L. Tretkoff, M.D. Tretkoff, Combinatorial group theory, Riemann surfaces and differential equations, 

Contemporary Mathematics, 33, 467-517 (1984). 
[37] M. Trott, Applying Groebner Basis to Three Problems in Geometry, Mathematica in Education and 

Research 6 (1): 15-28 (1997). 
[38] V. Vinnikov, Self-adjoint determinantal representations of real plane curves, Math. Ann. 296, 453-479 

(1993). 

[39] O. Ya. Viro, Gluing of plane real algebraic curves and constructions of curves of degrees 6 and 7, 
Topology (Leningrad, 1982), Lecture Notes in Math., vol. 1060, Springer, Berlin, pp. 187-200 (1984). 

Centre de recherches mathematiques Universite de Montreal, Case postale 6128, Mon- 
treal H3C 3J7, Canada 

E-mail address: kalla@crm.umontreal.ca 

Institut de Mathematiques de Bourgogne, Universite de Bourgogne, 9 avenue Alain 
Savary, 21078 Dijon Cedex, France 

E-mail address: Christian.Klein@u-bourgogne.fr 



